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Abstract 

We give an explicit formula for the solution to the initial value 
problem of the full symmetric Toda hierarchy. The formula is obtained 
by the orthogonalization procedure of Szego, and is also interpreted 
as a consequence of the QR factorization method of Symes ||. The 
sorting property of the dynamics is also proved for the case of a generic 
symmetric matrix in the sense described in the text, and generaliza- 
tions of tridiagonal formulae are given for the case of matrices with 
2M + 1 nonzero diagonals. 

Mathematics Subject Classifications (1991). 58F07, 34A05 



1 Introduction 

The finite non-periodic Toda lattice hierarchy can be written in the Lax form 



-^-L=[B n , L] , n = l,2,..- (1) 



1 



where L is a symmetric "tridiagonal" matrix, and B n is the skew symmetric 
matrix obtained from L by 

B n = l[L n = (L n ) >0 -(L") <0 , (2) 

a 

where (L n ) >0 , <q n denotes the strictly upper (lower) triangular part of L n . 

An immediate consequence of ([I]) is that the eigenvalues of L are invariant 
with respect to the time t = (ti, t2, • ■ •), and the Toda equation, correspond- 
ing to the case n = 1, as well as the hierarchy (|1|), can be shown to be a 
completely integrable system. Moser J| showed that for the Toda lattice 
(n = 1) the matrix L converges as t\ — > oo to a diagonal matrix contain- 
ing the eigenvalues of the matrix L(ti = 0), arranged in decreasing order. 
This sorting property has been further studied in O , where a general frame- 
work for calculating the eigenvalues of a symmetric matrix using the Toda 
equations was developed. 

In UJ, Deift et. al. extended the equation ([!]) to a system where L is an 
arbitrary symmetric matrix, which we call the full symmetric Toda equation. 
They then proved its complete integrability by constructing the constants of 
motion and the associated action - angle variables as a Hamiltonian system. 

In this paper we give an explicit formula for the solution to the initial 
value problem of the full symmetric Toda hierarchy, and show that the sorting 
property also holds in this case. The formulae we construct for the eigen- 
functions of L turn out to be an orthogonalization procedure described in 
Szego's book M. This result is also obtained via the QR factorization solu- 
tion method of |J. We also construct alternative formulae for the cases of 
symmetric 2M + 1 diagonal matrices, i.e. = for \i — j\ > M + 1. 

2 Full Symmetric Toda Hierarchy 

The full symmetric Toda hierarchy (|I|) can also be written as the compat- 
ibility condition of the set of two linear equations for the orthonormalized 
eigenvectors of the N x N symmetric matrix L, 



14 = \(j> , (3) 
d 



Let us here give some properties of the eigenvalue problem (|3|). Since the 
matrix L is symmetric, it can be diagonalized by an orthogonal matrix $ G 
0(N), i.e. 

L = $ diag ■■■ ,\ N )$ T . (5) 

Here the matrix $ contains the eigenvectors of L, say 0(A&) = (<f>i(\k), ■ • • , 
<MA fc )) T for fc=l,2,...,iV, 

$ = [^0, ••• , #A*)] = [^(A,)]^.^ • (6) 
Since $ G O(N), we also have the orthogonality relations 

JV 

0i(Afc)0j(A fc ) = % , (7) 

fc=i 

53 0fc(Ai)0fe(A i ) = % . (8) 
fc=i 

Then the entries of L, which we call be expressed as 

< A<fo</>j >= 51 A fc ^i(A fe )0j(A fc ) . (9) 



r/ 



fc=i 



Here the bracket <> is defined as a linear functional on F, the integrable 
functions of A, with the Dirac measure, i.e. 



< > : F -> R 

/(A)-/ R /(A)da(A) 



(10) 



where the measure da(X) is 



N 

da{X) = <K A - A fc )riA . (11) 
fc=i 

Note also that the entries of L n , which we will call a£j , are given by 

4° = (L n ) l3 =< X n Uj > • (12) 
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To solve the equations (0) and (|]), we first gauge transform the eigen- 
function 0(A,t) by the diagonal matrix 

A = diag( < >~ 1/2 , < ifj 2 N >" 1/2 ) , 

<p = Aip. (13) 

Note that the transformation ( |T3"D includes a freedom in the choice of ip, i.e. 
C [l3f) is invariant under ipi — > /^(t)^, with {/j}^ arbitrary functions of t. 
With dg), (0) and (§ become 

(A _1 LA)^ = , (14) 

A* = (A- W -(|-logA)^ (15) 

Writing 

(A _1 £ n A) = -2(A- 1 L n A) <Q + (A" 1 ™) -diag(L n ) , 
equation (|TJ) becomes 



"n 



-2 (A- X L" A) <o ^ + A> - (L n ) + A log A^ . (16) 



We here observe that ( |T6D can be split into two sets of equations by fixing 
the freedom in the determination of tp. In the components, these are 

iAlog<V?> = 4" 1 - (18) 
It is then easy to see that we have: 



Lemma 1 The solution of ftT7[ ) can be written in the form of separation of 
variables, 

^(A;t)=A(t)0°(A)e^ t ) , (19) 

where A(t) is a lower triangular matrix with diag(A(t)) = A(t = 0) = In, 
the N x N identity matrix, 0°(t) = 0(A;t) ; and f(A,t) = ££° =1 A n t n . 
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Note that we have chosen the initial conditions of if)(\, t) to coincide with 
0(A,t), ip(X;0) = 0°(A), and thus < ^ j > (t = 0) =< <f^] >= %. As a 
direct consequence of this lemma, and the orthogonality of the eigenvectors 
(0), i.e. < ipi^j >= for i ^ j, we have: 

Corollary 1 (Orthogonality): For each i G {2, ■ ■ ■ , N} , we have for all t 
with t m E R ; 



N 



< A<f>y {X ' t} >=E^(^;t)0°(A fc )e« A - t ) = 



(20) 



k=l 



for j = 1,2, • ••,«-!. 



Now we obtain the main theorem of the section, which yields the formula 
for the eigenvectors of L in terms of the initial data {</>°(A)}i<j<Ar: 

Theorem 2 T/ie quantities ^(A;t) are given 6?/ 

cu ... n, 

P«(A,t) 

^(A;t) 



A-i t 



#(A) 



• • Q— li 

•• $(A) 



(21) 



where c^(t) =< (fj^cj^e 2 ^^ >, and D k {t) is the determinant of the k x k 
matrix with entries Cy(t) ; i.e. 



D,(t) 



(c«(t)) : 



l<ij<fc 

Proof. From equation (^) with (|T9"D, we have 



(22) 



]>> jfe (t)<0 o fc 0°e 



>= , /or 1 < £ < i - 1 



(23) 



fc=i 



Solving (|23|) for A ik with = 1, we obtain 



Afe(t) 



A-i(t) 



(24) 
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where £)*_i(t) is the determinant A-i(t) in the form fl2"2|) with the replace- 
ment of the kth column (ci^, . . . , Cj_i fc) T by (cij, . . . , Cj_i j) T . From ([3~9|) , 
we have 



o 

h9k 



k=l 

3 e(A,t) y» 



+ e 5 



A-i 



e €(A,t) 

aT 



x -0? 



Cli C21 



Cli-1 



■ ■ + 



cn • • • Cu-l 



+ 



C]i-1 



which one can see is just Q21|) . □ 

We note from fl2~T|) that < ^ > can be expressed by A 



< >= 



A-i 



This yields the formulae for the normalized eigenf unctions 

Cn ... Ci 

P €(A,t) 

<MA,t 



A(t)A-i(t) 



Q— 11 ■ • • Cj— lj 
0?(A) ... $(A) 



(25) 



(26) 



This derivation of (|26| ) is nothing but the orthogonalization procedure 
of Szego, which is equivalent to Gram - Schmidt orthogonalization. Indeed, 
using the fact that fa = ,fii/2 , we see that equation ( |T9| ) expresses the 
relations 

fa(X,t)espan{^(\)e^ t \---A^)e^' t) } , < = 1, - - - , JV, (27) 

where we are viewing {4>i(\,t)}f =1 as a collection of orthogonal functions 
defined on the eigenvalues of the matrix L. Relations (p7[), together with the 
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orthogonality 

< fcfc >= Sy , (28) 

imply that {0«(A,t)}^ 1 are obtained by orthogonalization of the sequence 

|0°(A)e^ A ' t - ) |. , and hence one may obtain equation ( |2"E| ) from classical for- 
mulae, directly from (pTj) and fl2~8"|). 



With formula (p6|) , one can solve the initial value problem of (ff) and 

via 

Oy(t) =< A0i^ > (t) . (29) 

We also remark that the determinant Dk in (12) gives the tau-function 
of this hierarchy, and plays a fundamental role in integrable systems H. 



3 The QR Factorization Method 

Here we give an alternative derivation of the formula (^) by using the QR 
factorization method of Symes ||. 

Let L° be the initial matrix of L(t), i.e. L(0) = L°. Symes' factorization 
method is to compute the unique decomposition of the symmetric matrix 

e «Z°;t) = QF(t)R(t) (30) 

where f(L°;t) = En=it n (L°) n , Q(t) G O(N) with Q(0) = I N and R(t) is 
an upper triangular matrix with Ra > Then the evolution of L(t) can be 
obtained by the relation 

L(t) = Q(t)L°Q T (t). (31) 
The orthogonal matrix $(t) defined in (|6]) is then given by 

$(t) = Q(t)$° (32) 

where $° = (^(Aj))^*. 

Now from the diagonalization of L° in (|J), we have 

e e(L°;t) = ^(A i;t)) . . . ; e «A„;t)) ( $0 ) T (33) 
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which, together with ( |30| ) and ( |3~1"D leads to 



diag (e^\ • • • , e^^) ($°) T ^(t)" 1 = $ T (t) . (34) 



Equation ( |34|) should be compared with (|i9|); we note that ( |34J ) immediately 



implies relations (|27l). Hence, as described following relations (|27|) , the fac- 
torization method of Symes is enough to establish Theorem and in 
particular formula (p6|). 



4 Sorting Property 



Here we consider the case of the full symmetric Toda equation, that is, we 
have only one time t\, which we denote by t. As a consequence of the explicit 
formula for 0j(A; t) in Theorem ^ we can now show the sorting property, i.e. 
as t —>■ oo, the matrix L(t) converges to a diagonal matrix with ordered 
eigenvalues. In order to show this, we first show: 



Lemma 2 Let the eigenvalues of L be ordered as Ai > A2 > 
as t — > 00, we have 



D n (t)exp ^-2 htj -> (det$° n y 



> Xn- Then, 



(35) 



where $° is the n x n matrix given by $° 



Proof. From (|22|) , we have 



AT AT n 

^ = E---E n^)e 2t ^ 

Ai=l A n =l i=l 



0?(A! 



) (^i))l<j,j<n 

■•■ 0?(A„) 
■•■ <(A.) 

-JV 



D n (t) = E d{l)x^ 

\£\=n 



■ X 



N 



(36) 



Expanding this in terms of exp {2(£ } X)t} with (£, A) = J2i=i 4A« for ^ e Z> , 
and with \£\ = Y,i=i U — n i we obtain, setting Xj = exp (2Ajt), 

-4v 



(37) 
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We first show that the nonvanishing dominant term (for t large) in fl3"T|) is 
given by the case £ = (1, . . . , 1, 0, . . . , 0), that is, £{ = 1, for i = 1, . . . , n, and 
ii = otherwise. 

To show this, we note that 



d{£) = 



(38) 



if there exists £i > 2 for some i. Indeed, for such a term, the determi- 
nant in ([36|) would necessarily vanish. This, along with the ordering of 
the eigenvalues, implies that the dominant term in ( j37|) is the one with 
£ = (1, . . . , 1, 0, . . . , 0), and is given by 



rf(l,...,l,0,...,0) 



d" 



dx\ ■ ■ ■ dx, t 

Which is, summing over all permutations 



-D n (t) 



1 2 

£i £2 



n 



xi=---=xn=0 



(39) 



d(l,...,l,0,...,0) 



E 



0?(A £l )0?(A £l ), 
0i(A*J<(A,J, 

E0 1 (a. 1 )---<(a.„; 

Pn 



•• 0?(A, x )0°(A 4 ) 
0?(A,J ••• 0°(A £l 



<(A£ 



(40) 



Now we have 

Theorem 3 Suppose that det§^ 7^ for n = 1, . . . , N. Then as t — > 00, t/ie 
eigenf unctions <pi(Xj;t) satisfy 



(j)i(Xj) — > 5ij x sgn (det^fj sgn (^det^_^j 
which implies L(t) — > diag (Ai, . . . , Atv) ; /rom 



(41) 
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Proof. From the formula (^6|) and Lemma ^] we have, as t — > oo 
0n(A; t) 



|dei$° | |det$°_i| exp [(2 YJCi A; 



A n ) t 



X 



X 



Cll(t) 
Cln-l(t) 

0?(A) 



Cln(*) 



Cn— In 

<(A) 



The dominant term in the determinant gives 



0?(A) 



<(a^: 

0n(A) 



0?(Ax) • • • 

#(V-i) ... 

0?(A) • • • 

Pn-1 



0n(Al 



<(An-i; 

0n(A) 

••0°-i(A,„_ 1 )#({^}p 1 ) 



(42) 



(43) 



where #({£j}™ = l) is the signature of the permutation {tj}™ = \. We note that 

A,-, 3 = 1, • 



n-l 
J=l- 

n — 1. The ordering 



in — 1 
lj=l 

the determinant in (|4^) is zero for A 
Ai > ■ • • > A at implies the result. □ 

Observe that the arbitrariness in the selection of eigenvectors of the initial 
matrix appears in the asymptotic formula (f4T|). However, independent of this 
choice (naturally), we see from Theorem 3 that for any generic symmetric 
initial matrix, that is, Lq satisfies ^ for j = 1, . . . , N, and has distinct 
eigenvalues, the Toda lattice converges to a diagonal matrix, containing the 
eigenvalues arranged in decreasing order. 



5 Finite-Band Toda Hierarchy 

As a final comment, we give alternative formulae for 0j(A,t), for the case 
where the matrix L has (2M + 1) nonzero diagonals, M — 1 (tridiagonal) , 
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M = 2 (pentadiagonal), and so on. Thus we have a^- = for |z — j| > M. 
We further assume a^- 7^ for \i — j\ = M, and that the eigenvalues of L are 
distinct. 

For such a 2M + 1 - diagonal matrix, the first column of the eigenvector 
equation Lcj) = \cf) is 



M 



^aii^(A) +a x m+i0m+i(A) = A0i(A) . 



(44) 



i=i 



When the eigenvectors 0(A) are viewed as functions over the eigenvalues 
{Xj}jL 1: equation (|44]) expresses the relation 

0m+i(A) G span {^i(A), ... ,0m(A),A0i(A)} . (45) 

Likewise, we have 

0A/+i(A) G span{0i(A), . . . , </> M (A), A0i(A), . . . , A0j(A)} (46) 

for 1 < j < M. We leave the reader to verify that we have 

0fcM+j(A) G span {0i (A), . . . , 0m(A), A x [0i(A), . . . , 0m(A)] , . . . (47) 

. . . , A* x [0i (A), . . . , 0j(A)]j , 

A; = 0,1,... 
j = l,...,M • 

Relations (f47|), together with the orthogonality (p8|), yield the following or- 
thogonalization procedure. 

Define, for k = 0, 1, . . . and j = 1, . . . , M, 



and 



Then we have 



PkM+j = A fc 0j(A) 



Cij =< PiPj > 



(48) 
(49) 



&(A) 



Q-l 1 

Pi 



Pi 



(50) 
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where 



D n det [cij\ l<i - <n 



(51) 



We remark at this point that if ctk M+i , (k +i)M+i — for some 1 < i$ < 
M — 1, k > 1, then the procedure is modified by the replacement of pkM+i 
by \ k ~ ko (p ko M+i for k > k . 

Now from the preceding sections, when the matrix L(t) evolves according 
to the Toda hierarchy, we have the relations 



<Wi(A,t) e span{cf ) U\)e^ t \...,4M + / iKt} } 
for k > and 1 < j ' < M, which, together with fl47|) imply 



feM+ i(A, t) G span {A x [</>?(A)e^, . . . , < / (A)e«^ 



(52) 



(53) 



A fc x 



and hence for all t with tj > 0, we have the following orthogonalization 
procedure. Define 



p kM+j (\,t) = XWS™ 



fc = 0,... 
j = l,...,M 



and 



Then we have 



Cy(t) =< Pi/Oj > (t) 



cn(t) 



(54) 



(55) 



0i(A,t) 



Ci-l i(t) 
Pl (A,t) 



Cli(t) 
Cj-l i(t) 

/0i(A,t) 



(56) 



where 



D n = c/et [^-(t)]^^.^ 



(57) 
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